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We propose a new quantum model interpolating between the fully frustrated spin- 1/2 Ising model 
in a transverse field and a dimer model. This model contains a resonating-valence-bond phase, 
including a line with an exactly solvable ground state of the Rokhsar-Kivelson type. We discuss the 
phase diagrams of this model on the square and triangular (in terms of the dimer representation) 
lattices. 



I. INTRODUCTION 

The concept of the resonating-valence-bond (RVB) 
phasei has remained for many years one of the most 
fascinating topics in the physics of strongly corre- 
lated systems. This phase was originally proposed 
for frustrated spin systems^ and for high-temperature 
superconductors^ but to date no confirmed examples of 
physical frustrated magnets with the RVB phase have 
been found and the relevance of RVB physics to high- 
temperature superconductivity remains at the level of 
conjecture. However, many model systems have been 
proposed to exhibit the RVB phase^^"— including quan- 
tum dimer modelsf^i^ where this phase is most easily 
accessible. Interest in the model studies of the RVB 
phase has been further stimulated by experimentaliS and 
theoreticalii works on the spin-1/2 Heisenbcrg model 
on the kagome lattice, which suggested a spin-liquid 
phase. More recently, it has been suggested that this 
kagome-lattice model is closely related to a certain class 
of quantum dimer modelsJ^ In view of these develop- 
ments, model studies of phase diagrams of various sys- 
tems involving the RVB phase may be useful for further 
search and identification of this phase. In the present 
work, we propose a new class of models that realize an 
interpolation between quantum dimer models (QDMs) 
and fully frustrated spin-1/2 Ising models in a transverse 
field (FFIMs) with an exactly solvable RVB ground state. 

There is a well-known correspondence between QDMs 
and FFIMsi ^^i^^ When these two models are formulated 
on dual lattices (so that the sites of the Ising lattice cor- 
respond to the plaquettes of the dimer lattice and vice 
versa), frustrated and unfrustrated bonds of the Ising 
model can be put into correspondence with the presence 
and absence of a dimer, respectively. Therefore, in the 
limit of a strong Ising coupling (compared to the trans- 
verse field) , the number of frustrated bonds at each Ising 
plaquette should be odd and minimal, which implies the 
QDM constraint of fully packed dimers. 

An exact mapping between the two models exists, how- 
ever, only in the above-mentioned limit of the strong Ising 
coupling. In the well-studied examples of square and 
triangular latices, this limit belongs to the crystallized 
phase (dimers or Ising spins order and break the transla- 



tional symmetry of the lattice) On the other hand, 

in the opposite limit of a strong transverse field FFIM is 
in the disordered phase quite similar to the RVB phase 
realizable in QDMs on nonbipartite lattices. The char- 
acteristic features of the RVB phase are unbroken trans- 
lational invariance (on lattices with a single site per unit 
cell), exponential decay of correlations, topological order 
on multiply connected domains, and Z2 vortices (visons) 
as elementary cxcitationsi^^r— 

Although a conjecture was made that the disordered 
phase of the FFIM is continuously connected to the RVB 
phase of the QDMj^i^ this never was demonstrated ex- 
plicitly. In the present work we provide a justification 
for such an identification of the disordered FFIM phase 
and the RVB dimer phase by constructing a more gen- 
eral model which realizes a continuous interpolation be- 
tween FFIMs and QDMs. This model is exactly solvable 
(in terms of its ground state) along a line belonging to 
the RVB phase and connecting the high-field limit of the 
FFIM to the Rokhsar-Kivelson point of the QDM. 

The paper is organized as follows. In Sec. HIl wc briefly 
review the QDM and the FFIM and their RVB phases. In 
Sec. mil we construct the interpolating model connect- 
ing the QDM and the FFIM on mutually dual lattices. 
In Sees. IIVI and fVl we analyze in more detail the phase 
diagrams of this interpolating model on the square and 
triangular (in terms of the dimer representation) lattices, 
respectively. Finally, in Sec. IVIl we summarize and dis- 
cuss our findings. 



II. RVB STATES IN DIMER AND SPIN 
MODELS 

A. Quantum dimer model 

The simplest QDM is the so-called Rokhsar-Kivelson 
(RK) dimer modeli^ 
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where the sum is taken over all tetragonal plaquettes of 
the lattice. At v = t, both on square and triangular lat- 
tices, the ground state of the RK model (or, possibly, one 
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of the ground states, depending on the topology of the 
lattice cluster) is given exactly by the sum of all possi- 
ble dimer configurations with equal amplitudes (the RK 
state) ,^ 

On the triangular lattice, the RK state is an exem- 
plary RVB state: it has a finite gap and exponentially de- 
caying correlations, topological degeneracy and visons as 
excitationsi^^"— It however lies at the phase boundary: 
at V > t, the model crystallizes via a first-order phase 
transition into the so-called staggered (or nonflippable) 
state. On the other side, at v < t, there is a finite win- 
dow of RVB states. According to numerical studies,—!^ 
this window extends to approximately v/t = 0.83(2). 

On the square lattice, a,t v — t, the RK state is critical: 
it has gapless excitations and power-law correlations. On 
both sides of the point v = t, the system crystallizes as 
soon as one deviates from this point J^i^Sr— 

The difference between the square and triangular lat- 
tices results from one of them being bipartite while the 
other is not. There is a vast literature on the phase di- 
agram of the RK dimcr models on both lattices and on 
the properties of the RVB states (for a review, see Rcf. 

B. Fully frustrated Ising model 



The Hamiltonian of the FFIM reads 
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where the coefficients Mij = ±1 are chosen in such a way 
that their product over any plaquette of the Ising lattice 
is negative (with both parameters J and T assumed to 
be positive) Here and below we denote the positions 
of Ising spins by Latin indices i, j, etc., whereas the first 
sum in Eq. Q is taken over all pairs of nearest neighbors 

In this model, the RVB state is obviously realized in 
the limit of strong transverse magnetic field T ^ J. In 
this limit, the ground state has all spins almost fully po- 
larized along the field, their transverse components being 
disordered. Since each spin aligned in the x direction is 
a linear combination of spins cr^ = ±1 with equal am- 
plitudes, the resulting ground state contains all sets of z 
projections of spins with almost equal amplitudes;^ which 
resembles the RVB state in the RK model. This state is 
translationally invariant, has only local correlations, and 
has an excitation gap of order F. The Ising spins in the 
FFIM model correspond to the vison operators in the 

C. Relation between the two models 

A rigorous mapping between the FFIM and the QDM 
exists only in the zero-field limit (F/J 0) of the 
FFIM when it becomes equivalent to the RK model with 
v/t = OJ^ This mapping involves FFIMs and QDMs on 



mutually dual lattices (i.e., plaqucttes of one lattice cor- 
respond to sites of the other) . While the construction of 
the interpolating model developed in the next section is 
generally applicable to any dual pair of lattices, we will 
further illustrate it with specific examples of the triangu- 
lar and square QDM lattices (which corresponds to the 
hexagonal and square FFIM lattices, respectively). On 
both of these lattices, the point at which there exists a 
rigorous mapping between the FFIM and the QDM cor- 
responds to crystallized phasesJ^"— 

Below we explicitly construct a model which realizes 
a continuous interpolation between the FFIM and the 
QDM on mutually dual lattices and is exactly solvable 
along the line connecting the FFIM at F/J oo (the 
limit of noninteracting spins) to the RK model at v = t. 
Along this line the system belongs to the RVB phase. 



III. INTERPOLATING MODEL 

A. Construction of an RK-type model 

For the interpolating model, we use the same Hilbert 
space as in the FFIM. It will be convenient to introduce 
the basis Z defined in terms of projections of spins on 
axis z. Then we first postulate the ground state of the 
interpolating model parametrized by a "chemical poten- 
tial" fi: 
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where the sum is taken over all Ising-spin configurations 
|c) from the basis Z. To avoid confusion, the lattice at 
whose sites the variables aj are defined is called below the 
FFIM lattice and the lattice dual to it the QDM lattice. 
The sites of the latter are denoted by Greek letters. Note 
that the construction described in this section is rather 
general and does not require the lattices to be periodic. 
In terms of the dimer representation, the quantity 



(4) 



describes the number of dimers Uap = (1 — Tij)/ 2 on the 
bond a/3 of the QDM lattice which crosses the bond ij 
of the FFIM lattice (i.e., Tij equal to plus or minus one 
corresponds to the absence or presence of a dimer, re- 
spectively). Because of the frustration of the parameters 
Mij, each site of the dimcr model lattice must belong to 
an odd number of dimers. In the limit /i — > cx) only con- 
figurations with exactly one dimcr per site survive, thus 
the ground state ^ continuously interpolates between 
the state which is fully polarized in the x direction (the 
ground state of FFIM in the limit T / J ^ oo) at /i = 
and the RK state (the ground state of RK dimcr model 
&t V = t) aX ^ ^ oo. 

At the next step, we construct the Hamiltonian whose 
ground state is given by ([3]). This is performed in the 
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"supersymmetric" way suggested by Henley^i; namely, 
the Hamiltonian is assumed to have a form of a quadratie 



H 



(5) 



where the operators Qa are such that they annihilate the 
chosen ground state. We call such a class of models the 
RK-type models. 

The decomposition ([5]) for the standard RK model 
([T]) is well known, with the operators Qa removing two 
dimers on one tetragonal plaquette a in two possible ways 
with opposite signsi^ In a shorthand notation, this op- 
erator may be written as 
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Note that this operator does not acts in the physical 
Hilbcrt space but maps physical configurations of dimers 
onto configurations in some auxiliary space with one pla- 
quette removed. 

The above construction can be directly generalized to 
Ising-type systems to annihilate the ground state 
The simplest choice of operators Qa involves removing 
one Ising spin. 



Q^ 
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where the sum over j is taken over the nearest neighbors 
of the site i denoted n.n.(i). 

After some simple algebra, the Hamiltonian ([5]) with 
operators Qi given by ([7]) can be rewritten as 



exp 
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where the sums over i are taken over all sites of the FFIM 
lattice. The sum over j for each i is taken over the nearest 
neighbors of i. 

In terms of the dimer representation the first term in 
Eq. ([8|) can be rewritten as J^i "^xp [—fJ.{mi — 2ni)]. Here 
nii is the number of nearest neighbors of site i on the 
FFIM lattice (i.e., the number of bonds belonging to 
Vi, the corresponding plaquette of the QDM lattice) and 
Hi = X]{a^}6'P "'a/5 total numbcr of dimers on Vi- 

In the case of a simple periodic lattice rrii = const. At 
the same time, operator af corresponds to the inversion 
of dimer occupation numbers (uap i— )■ 1 — Uap) on all 
bonds belonging to Vi. 

The resulting Hamiltonian is an RK-type model 
which continuously interpolates between the system of 
non-interacting spins in a uniform magnetic field [in other 
terms, the F/J ^ cx) limit of the FFIM ([2])] at ^ 
and a dimer model with the RK ground state at /i — > oo. 
More precisely, in the /i oo limit, the model ([5]) splits 
into sectors with different number of dimers (including 



possible overlaps of dimers), and it is the sector with 
the minimal number of dimers (i.e., with non-overlapping 
dimers) which contains the RK ground state 

Furthermore, we can show that on the whole line 
< ^ < oo the ground state ^ has a finite correla- 
tion length (at least, on the commonly used lattices). It 
follows from the observation that the equal-time correla- 
tions in the ground state ([3]) are exactly the same as in 
the classical fully frustrated Ising model without the field 
defined on the same lattice at the dimensionless temper- 
ature T = For quite a number of two-dimensional 
lattices the properties of such models are known from the 
exact solutions (for a review, see, e.g., Ref. [H). In par- 
ticular, the models on square and triangular lattices are 
critical at T = and acquire a finite correlation radius 
at an arbitrarily low temperature, whereas the models on 
honeycomb, kagome and pentagon lattices have a finite 
correlation radius already at zero temperature. In all 
these cases the correlation radius continuously decreases 
with the increase in T and at T — > oo shrinks to zero. 

We can also claim that this model (at < ^ < oo) has 
topological order. Indeed, different topological sectors 
of the dimer representation in terms of the Ising rep- 
resentation correspond to different boundary conditions 
for spinS)^ and the absence of translational symmetry 
breaking in combination with a finite correlation radius 
guarantees that these sectors are degenerate in the ther- 
modynamic limit, as expected in a topologically ordered 
state. 

Note that for n < oo this topological order is of the 
Z2 type regardless of the lattice geometry (e.g., also for 
the square lattice) . This contrasts with the properties of 
the pure RK dimer model, which is of the U(l) type on 
bipartite lattices, i.e., incorporates an infinite number of 
topological sectors which can be characterized by inte- 
ger "winding numbers." The reason for this reduction of 
symmetry is that the U(l) conservation law for the dimer 
model on the square lattice breaks down to Z2 as soon 
as non-dimer states are allowed. A Z2 topological or- 
der implies also the existence of vortexlike Z2 excitations 
(visons).— In fact, vison excitations in the Ising-spin lan- 
guage are generated by the erf operatorsi^i^ 

While the ground state of our model is exactly 
known for any lattice, the excitations are not. How- 
ever, since the model is of the RK type, then, as 
pointed out by Henley, the spectrum of the lowest ex- 
citations can be efficiently computed with the classical 
Monte Carlo methodf^i^ In principle, the spectra of 
both vison and non-vison excitations can be computed 
by modeling a classical stochastic walk in the space of 
configurations 



B. Generalized two-parameter model and its 
reduction to the FFIM 

A continuous interpolation between the FFIM with an 
arbitrary ratio F/J > and a QDM with a variable pa- 
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rameter describing the relative strength of different terms 
can be achieved by a shght modification of the Hamilto- 
nian ([8]). It consists in ascribing independent amphtudes 
to the potential and kinetic terms, 



H{U,W) = C/^cxp 



j"=n.n.(i) 



Obviously, only the dimensionlcss ratio of U and W is of 
relevance for the phase diagram. 

This model reduces to an FFIM in the limit when fi 
is taken to at a constant value of the product J7/i. In 
this limit, one can expand the exponent in the first term 
of Eq. (|9]) and keep only the terms linear in /i, because 
the coefficients in front of all higher order terms vanish. 
The summation of contributions from different plaquc- 
ttes then immediately leads to the FFIM Hamiltonian 
([2]) with J = and T = W. This approach is univer- 
sal in the sense that it works for any lattice. 

The reductions of the model © to dimer models are 
more delicate and have to be discussed separately for 
different types of lattices. Below wc explicitly analyze 
the model (jH]) on the square (Sec. IIV|) and honeycomb 
(Sec. [V]) lattices, which correspond to dimer models on 
square and triangular lattices, respectively. 



IV. SQUARE LATTICE 

On the square lattice, quantum dimer models can be 
obtained from the model ^ in two different limits. 



A. QDM at ^ ^ 

One possible reduction to the standard RK model ([T]) 
is obtained from ^ by taking the limit /.t — > in a way 
different from that described in Sec. IIII Bl At < /i <C 1, 
the dimer states (in which each site belongs to only one 
dimer) are separated from all other states (not satisfying 
this rule) by the gap of the order of U fj,. Therefore, if U 
goes to infinity faster than fi goes to zero, the energies 
of all states with overlapping dimers go to infinity and 
the Hilbert space of the system is reduced to that of 
the dimer model. In particular, when /x — > at U^^ = 
const, only the second-order terms in the expansion of 
the exponential remain finite. One can show that their 
contribution to the potential energies of different dimer 
states is proportional to the number of square plaquettes 
populated by two parallel dimers, which means that the 
potential energy acquires the form of the first term in 
Eq. ((J) with V = AUfi'^. 

When only the dimer states are allowed, the action of 
the kinetic term from Eq. ^ on them is reduced to the 
possibility of dimer flipping [like in Eq. ([T])] with ampli- 
tude t = W. Thus, when fi tends to zero at a finite value 
of U^^, the model © is reduced to the RK model ([T]) 
with V = 4U^i'^ and t = W. 



B. QDM at ^ ^ c» 

Quite remarkably, another reduction to a dimer model 
is implemented in the opposite limit fi ^ oo (at finite val- 
ues of U and W). In this limit, potential energy vanishes 
for all plaquettes without dimers and for plaquettes with 
one dimer, is equal to U for plaquettes with two dimers 
and is infinite for plaquettes with three or more dimers. 
This splits the Hilbert space of the system into sectors 
with different number of dimers, because all processes 
changing this number are prohibited. In particular, the 
sector with the smallest possible number of dimers corre- 
sponds to the standard situation when any site belongs 
to only one dimer. Within this sector the Hamiltonian 
(|9]) is completely equivalent to the RK Hamiltonian ([T]) 
with V = U and t = W. 

Other sectors involve configurations with overlapping 
dimers (see Fig.[T]) that we call below "stars." Since each 
star adds one extra dimer, these sectors can be labeled by 
the total number of stars in the system. We can therefore 
refer to the model obtained in the — )■ cxj limit as a "star- 
dimcr model." In terms of the height representation for 
dimer coverings of the square lattice^^i^ each star cor- 
responds to a screw dislocation, whose sign depends on 
to which of the two sublattices this star belongs. 

A relevant question in such a situation is whether the 
"starless" sector is indeed the lowest-energy sector of the 
model. At U = W, one can show that, while a single- 
star configuration has the same zero energy as the RK 
ground state, starting from two stars, the ground-state 
energy becomes positive. This is a consequence of a con- 
tact repulsive interaction between stars: when stars are 
close to each other, some flips on plaquettes with two 
dimers would produce plaquettes with three dimers and 
therefore are prohibited. Accordingly, the kinetic energy 
cannot fully compensate the potential one. Since on a bi- 
partite lattice stars are "charged" (with each star occupy- 
ing three sites from one sublattice and only one from the 
other), they can appear in a finite system only in pairs, 
and therefore the one-star sector can be disregarded. We 
can then summarize our conclusion that, at U = W , the 
starless sector has the lowest ground-state energy. 

Is is also possible to prove that the starless sector gives 
the lowest-energy state a,iU > W. In this case, the stag- 
gered state, in which each plaquette contains only one 
dimer and therefore no fluctuations are possible, real- 
izes the ground state of the model [the same nonflippable 
state is known^i to be the ground state of the RK model 
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FIG. 1: (Color online) "Star" configurations of dimers on the 
square lattice and their possible hopping. 
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FIG. 2: (Color-online) A schematic representation of the 
phase diagram of the two-parameter model ((9} on the square 
lattice, in the coordinates 



([T]) a.i V > t\. Indeed, since the full Hamiltonian may be 
represented as 

H{U,W) = H{W,W) + H{U -W,Qi), (10) 

and the nonflippablc state is a ground state of both 
H{W,W) and H{U - W,Q) (at U > W), it must also 
be a ground state of H{U, W) (without regard to a sec- 
tor). 

For U < W, the situation is more complicated: the 
sectors with stars may, in principle, have lower energy 
than the ground state of the RK model. While the lat- 
ter is known to be crystal-ordered (with the plaquette 
or columnar type of order— li^i^^i^) , it may be possible 
that in some interval of Y stars are energetically favor- 
able (which would possibly modify or destroy the crystal 
order). A simple perturbative study of stars doped into 
the liquid at the RK point Y = 1 suggests that stars are 
not energetically favorable in the star-dimer model near 
the RK point. However, for a rigorous justification of 
the absence of stars, a more careful (possibly numerical) 
study is necessary, which takes into account the presence 
of a crystalline order. We do not address this issue in the 
present work, but leave it for future studies. 



C. Phase diagram 

Constructing a full phase diagram of the model (|9]) 
is a challenging problem. The model depends on two 
dimensionless parameters /x and U/W, and already the 
boundary of this parameter domain contains several in- 
teresting models, including those discussed above (FFIM, 
RK dimer model, and the star-dimer model). 

For a better visualization of the discussed limiting 
cases, we choose the coordinates for the phase diagram 



U 
W 



sinh2(2M), Y 



W 

U 



coth(2^). (11) 



Our proposal for the phase diagram in these coordinates 
is schematically shown in Fig. [2l The vertical axis X = Q 
corresponds to the FFIM with T/J = Y, and the hor- 
izontal axis y = to the RK model with v/t = X. 
The X — > oo limit corresponds to the star-dimer model 
whose starless sector is given by the RK dimer model 
with t/v — Y. At least for < F < 1, the true ground 
state of the model ^ belongs to this sector. 

In terms of the original variables U, W and fi, the 
point X = oo, Y = corresponds to = and, natu- 
rally, can be achieved at any ^ (which in terms of Fig. [2] 
corresponds to approaching this point from any direction 
within the square). For any /i the potential energy of the 
considered model is minimized in the staggered states in 
which all plaquettes contain exactly one dimer, = 1. 
Comparison of this energy with the potential energy of 
the states with small numbers of plaquettes with Ui ^ 1 
suggests that for W/U <C max(l, /i) the quantum fluctu- 
ations on the background of a staggered state are weak 
and cannot destroy it. In conjunction with knowing that 
the staggered states are the ground states of the model 
both at X > 1, Y = and at X = oo, Y < 1, this allows 
one to expect that in a finite region of the phase diagram 
the ground state has the same type of ordering. 

The "RK line" U = W at which the ground state is 
given by Eq. ^ corresponds to F = 1 + 1/ \/X. Since at 
this line the system is in the RVB state with a finite cor- 
relation radius, one expects that the RVB state occupies 
a finite region on the phase diagram (in the vicinity of the 
RK line) . This region has to extend to the upper part of 
the axis X ^ 0, which corresponds to the translationally 
invariant disordered phase of the FFIM. 

Below a certain critical value of Y, the FFIM is in an 
ordered phase with a broken translational invariance. At 
r/ J = (which corresponds to the point X = y = in 
the lower left corner of our phase diagram), this phase is 
connected to the crystal phase of the RK dimer model at 
small values of v/t. 

Furthermore, one can argue that the crystalline orders 
present in the RK dimer model (i.e., at F = 0) extend 
to a finite region of the phase diagram at F > 0. In- 
deed, at small Y, the potential energy of a star is large, 
Est « iUiJ. K 2t/Y. Therefore, at F < 1, star config- 
urations appear only as virtual pairs and may be taken 
into account as perturbative corrections within the quan- 
tum dimer model. In particular, the lowest order pro- 
cesses lead to the rcnormalization of v and to the ap- 
pearance of an additional kinetic term related to a cyclic 
flip of three dimers belonging to two neighboring plaque- 
ttes. The amplitude of these corrections is proportional 
to 1/F, and therefore their presence cannot destroy the 
crystalline phases of the RK model in some region above 
the F = fine. 

As revealed by numerical studies of the RK dimer 
modelfi^i^^iSl at < < 1 its phase diagram contains 
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at least two different crystal phases: the plaquette and 
columnar ones. Therefore, the lower left part of our phase 
diagram must also contain two or more crystal phases. 

It follows from the analysis of Ref. [l^ that the transi- 
tion between the rightmost of these phases and the stag- 
gered crystal occupying the lower right part of the phase 
diagram may occur along one of the two scenarios: ei- 
ther as a first-order transition line terminating at the RK 
point (X = 1, y = 0) or as a devil's staircase (complete 
or incomplete) of intermediate commensurate phases. In 
either case, we do not expect any RVB region in the 
vicinity of the RK point. 

Yet another crystal region is expected to exist at the 
vertical axis X = oo above the RK point Y = 1. As men- 
tioned in the previous subsection, a simple variational 
analysis of stars doped into the RK state indicates that 
they should be energetically unfavorable just above the 
Y ~ 1 point. Therefore a plaquette crystal is expected 
in the dimer-star model adjacent to this point. 

A deviation from the X = oo axis leads to the mixing 
between the different sectors of the star-dimer model. 
The difference from the region just above the line Y ^ 
is that, in the vicinity of the X = oo axis, the proper 
energy of a star is not high. However, at X ^ 1, the 
amplitude of the formation of a pair of stars is low. In 
addition to that, the presence of the crystalline order in 
the starless sector induces a linear in distance attraction 
between the stars. It appears because, in terms of the 
height representation, stars correspond to screw dislo- 
cations. The combination of these factors leads to the 
confinement of stars and does not allow them to destroy 
a crystalline order in some vicinity of the X = oo axis 
(more precisely, of its part where the star-dimer model is 
in the crystal phase). On the other hand, when deviation 
from the X = oo axis takes place at y = 1 , the confining 
interaction between the stars induced by the crystalline 
order is absent and one immediately gets into the RVB 
phase. 

Further interesting phases may be possible in the upper 
right corner of the phase diagram, with stars playing a 
role in the energetic balance. In the present work, wc do 
not explore this region of the phase diagram, but leave 
this for future studies. 

V. TRIANGULAR QDM / HONEYCOMB FFIM 
LATTICE 

In this section, we consider the same model © in the 
case of the triangular QDM lattice, which corresponds 
to the honeycomb FFIM lattice. As on the square lat- 
tice, the FFIM limit is achieved in the limit ji ^ 
taken at Ufi/W ~ const, as is explained in more de- 
tail in Sec. IIII Bl The dimcr limit is, however, rather 
different from the square-lattice case. The reason for 
this difference is that, on the triangular lattice, the con- 
straint of nonoverlapping dimers prohibits direct flips in- 
duced by operators cr| , and the dynamics of dimers has to 
be mediated by virtual fiips via intermediate non-dimcr 



states. Since the magnitude of the gap to the nondimer 
states in the Hamiltonian © is of the order of ?7sinh/z, 
the dimer model can be expected to be realized when 
PF/(?7 sinh/i) — ^ 0. In this limit, the only terms which 
survive in the effective Hamiltonian for the dimer model 
come from the second order of the perturbation theory. 

Similarly to the "star-dimer" limit on the square lat- 
tice, the peculiarity of the limit W/ [U sinh /i) — )• on the 
triangular QDM / honeycomb FFIM lattice is that the 
potential term in Hamiltonian ^ does not impose the 
rigorous dimer constraint but prohibits only having more 
than one dimer on every triangular plaquette. In addi- 
tion to the usual (non-overlapping) dimer coverings this 
also permits the star-like overlaps of three dimers form- 
ing angles 120° with respect to each other (sec Fig. 
As a consequence, the effective model obtained in the 
considered limit contains not only the dimcr sector, but 
also the sectors with such three-dimer stars. As in the 
star-dimer model on a square lattice these sectors differ 
by the number of dimers: each star adds one extra dimer. 

A. Rokhsar— Kivelson limit 

We first discuss the RK-point limit, which is realized 
when W/{U sinh/i) — >■ while /x — >■ oo. In this limit, the 
only terms which survive in the effective Hamiltonian for 
the dimcr model arc related to virtual fiips on triangular 
plaqucttcs with one dimer. Analogous fiips on plaquettes 
without dimers can be neglected since they lead to inter- 
mediate states with three dimers per plaquette (which 
for fi ^ oo are infinitely higher in energy than the in- 
termediate states with two dimers per plaquette). Then, 
in the dimer sector (without stars), we arrive at the RK 
Hamiltonian with v = t = / {2U sinh ^) . 

In the case of sectors with stars, however, second-order 
perturbation theory shows an imbalance between the po- 
tential and the kinetic terms. Namely, for each star, 
there is a positive potential contribution of the order 
W'^/{U sinh/i), which is not compensated by a kinetic 
term. We then conclude that configurations with stars 
are separated by a finite gap from the RK dimer sec- 
tor, the magnitude of the gap being of the order of the 
coupling constant of the RK model. 

B. QDM at finite fj, 

A more general QDM is obtained when the limit 
W/{U sinh n) is taken at a finite value of fi. In this 




FIG. 3: (Color online) "Star" configurations of dimers on the 
triangular lattice and their possible hopping. 
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fo,o= [e2|^+ 1 +4e-2H]"'' 



fo,i = [e2^i + 2 + 3e-2|J] 




> fo,2=[e2n + 3 + 2e-2ti]"'' 




fo,3=[e2ji + 4 + e-2H]" 



fl,o=[1 +e-2ji]- 



fi 1 = 1/2 



FIG. 4: (Color online) The coefficients fn,fi entering Eq. (fT2|) 
for different configurations of dimers on the plaquette Vi 
(shaded triangle in the middle) and the three neighboring 
plaquettes. 



case, the effective kinetic term is also of the RK type 
with t = W^/{2Usmh^). The effective potential energy 
becomes, however, more complicated, because at < oo 
the virtual flips on triangular plaquettes without dimers 
also have to be taken into account. 

In any configuration with only nonoverlapping dimers, 
a virtual flip on the plaquette Vi produces a gain of the 
potential energy 



E,=-- 



2U sinh/i 



frii.rii (m) 



(12) 



where factors (m) depend on rii, the number of 

dimers on Vi, and fii, the total number of dimers be- 
longing to the three neighboring plaquettes of Vi but not 
belonging to Vi. Fig. U] shows the values of these fac- 
tors for all possible configurations of dimers on Vi and 
the three neighboring plaquettes. It is not hard to notice 
that 1/ fn,n{tJ-) depends linearly on both n and which 
allows us to replace the six formulas shown in Fig. Hlby 
a single expression. 



(13) 



be alternatively described in terms of local interactions 
between dimers. In this description, in addition to v, 
the energy of the interaction of nearest-neighbor dimers 
[see Fig. [SJa)], one needs to introduce two other cou- 
pling constants describing the interaction of next- nearest- 
neighbor dimers. We denote the pairwise interaction of 
such dimers as V2 [see Figs. [5l[b) and[5l[c)], whereas K 
denotes an additional three-body interaction in a loop of 
next-nearest-neighbor dimers [Fig. [5][d)]. These defini- 
tions imply that the dimer configurations shown in Figs. 
\Mjo) and [5jc) are ascribed the energy V2, whereas the 
configuration of Fig. [SJd) is ascribed the energy 3v2 + K. 
After straightforward algebra one obtains 

v/t = tanh^-2(/o,o-/o,i)> (14a) 
V2/t = -(/o,o-2/o,i + /o,2), (14b) 
K/t = /o,o - 3/o,i + 3/o,2 - /o,3 , (14c) 

where the argument of the functions /o,ft(/i) is omitted. 

In Eq. (|14ap . the relative magnitude of the second (neg- 
ative) term never exceeds 0.12 and vanishes in the limit 

— > oo. The functional dependence of v/t on ii is mono- 
tonic, with v/t ranging from to 1 as /i ranges from zero 
to infinity. Since the magnitudes of coupling constants V2 
and K are always numerically small (—0.0042 < V2/t < 
and < K/t < 0.0012, for any /z between zero and in- 
finity, see Fig. El), the resulting QDM can be understood 
as a small deformation of the RK Hamiltonian with v/t 
given by Eq. (jl4a[) . At ^ — and at /i — > oo both V2 
and K tend to zero and one recovers the RK model ^ 
with V = and v ~ t respectively. 

As for the star configurations, we have not analyzed 
their energies between the two limits ^ — > and ^ — )■ oo. 
Since at /i — >■ oo their energy is positive and at /i — > it 
becomes infinite, we conjecture that it remains positive 
also at intermediate values of /i, so that the dimer sector 
is always the lowest-energy sector of the model. How- 
ever, a rigorous verification of this property is probably 
possible only numerically. 



C. Phase diagram 

Based on the above consideration of the two limits, it 
is convenient to plot the phase diagram of the model (j9]) 




The potential energy given by the sum Ei may 



FIG. 5: (Color online) The potential energy given by the sum 
Ei may be reproduced by attributing energy i; to each 
pair of neighboring dimers (a), energy 112 to each pair of next- 
nearest-neighbor dimers [(b) and (c)] and an additional energy 
K to each three-dimer loop (d). 
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K/t 




FIG. 6: V2 and K as functions of v. 



on the triangular lattice in the coordinates ^ and 

W 

see Fig. [71 The vertical axis (/i ~ 0) corresponds to the 
FFIM with r/J = y, and the horizontal axis (y = 0) 
to the QDM. This QDM differs slightly from the RK 
model with v/t k, tanh ^ due to the presence of the ad- 
ditional interactions between the next-to-nearest dimers 
described by the coupling constants V2 and K [see Eqs. 
(|14bp and (|14cp ]. However these interactions are always 
very weak and vanish when /i — > and /i — > oo when the 
model ^ becomes (in its dimer sector) exactly equiva- 
lent to the RK model ([T]) with w = and v = t respec- 
tively. 

At v/t = 0.83(2) the RK dimer model on the tri- 
angular lattice experiences a phase transition from the 
RVB to the crystal phased A possible structure of this 
crystal phase has been discussed not only in the con- 
text of the RK dimer model but also in the context of 
the FFIM on the honeycomb lattice. The most likely 
crystal configuration has a 12-site unit cell in the QDM 
formulatio n^^i^^ corresponding to a 24-site unit cell in 
the FFIM formulation)^ We expect that, similarly to 
the RK model, our limiting dimer model obtained in the 
limit 2/ — > behaves in the same way. It seems to us un- 
likely that the very weak additional interactions by which 
this model differs from the RK model can lead to the for- 
mation of other (more complicated) crystals. Since the 
crystal phase with the same structure exists also in the 
FFIM at small enough values of y, it can be expected to 
occupy a finite region in the lower left corner of the phase 
diagram (see Fig. [7]). 

The RVB phase of the RK model spans the interval 
between v/t w 0.83 and v/t = 1, which corresponds 
to /i > 1.22, according to Eq. (|14ap . In this region, 
the next-nearest-neighbor interactions are very weak. In 
particular, at the transition point —V2/t ~ 10"'^ and 
K/t ^ 10^^ with a further rapid decay at v/t — > 1, the 
asymptotic behavior being V2/t k, — (1/4)(1 — v/t)^ and 
K/t ~ (3/8)(l — v/t)'^. This suggests that the presence 
of the additional interactions cannot induce a noticeable 
shift of the position of phase transition in comparison 



with the RK model or have any influence on the proper- 
ties of the RVB phase. 

The exactly solvable line is given in our coordinates 
by y = 1/(2 sinh fi) . It follows from the exact solution of 
the classical antiferromagnetic Ising model on the trian- 
gular lattice^ that the ground state ^ has a finite cor- 
relation length everywhere on this line, and therefore we 
expect that this line lies fully outside the crystal phase 
in our phase diagram. Note that the staggered (non- 
fiippable) phase of the RK dimer model, which appears 
at u/i > 1, is not present in our phase diagram, since 
this region of parameters cannot be reached within our 
model. We can only span the region between v/t = 
and v/t ~ 1, with additional interactions of next- nearest- 
neighbor dimers which are always too weak to stabilize 
the staggered phase. 

VI. CONCLUSION 

In this work, we have constructed a model interpolat- 
ing between a quantum dimer model and a fully frus- 
trated Ising model in a transverse field and studied its 
phase diagram for square and triangular (in terms of the 
dimer representation) lattices. This interpolating model 
exhibits an RVB phase, including a special "RK line" 
with an exactly solvable ground state. This RK-type 
ground state generalizes the construction at the RK point 
of the quantum dimer model, with its equal-time corre- 
lations given by the classical fully frustrated Ising model 
(without a field) at a finite temperature. 

While the RK dimer model is a well-known "proto- 
type model" of the RVB phase, our construction may be 
considered as its generalization including "star" configu- 
rations. It can be useful for studying the role of spin and 
charge degrees of freedom in RVB states. 




FIG. 7: (Golor online) A schematic representation of the 
phase diagram of the two-parameter model ((9]) on the tri- 
angular lattice in the coordinates fi and y = W/{2U sinh/i). 
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Indeed, a qualitative description of the RVB state 
in a spin system contains elementary spin excitations 
(spinons) and vortexlike excitations (visons). In ap- 
plication to high-temperature superconductivity, it has 
been also suggested that doping RVB liquids gives rise 
to charge excitations (vacancies or holons);^^""— At the 
same time, the RVB state in the pure dimer model only 
contains visons as elementary excitations, with spinons 
and holons prohibited by the close-packing dimer con- 
straint. From this point of view, the model constructed 
in the present work may serve as a prototype model of an 
RVB state including both visons and holons (since star 
configurations carry extra charge). 

Finally, we would like to note a certain similarity be- 
tween the appearance of RVB states on the square lattice 
close to the critical RK points in the model of the present 
work (the point X = od, F = 1 in the phase diagram in 
Fig. (2) and in that of Ref. [s^- In both situations, the 
initial RK state is deformed into an RVB state by intro- 
ducing a small concentration of charged objects (stars in 



the present work and diagonal dimers in Ref. IbsI ). One 
may note that, in both cases, the charged defects cost 
initially zero energy and the deformation is achieved by 
tuning the amplitude of their creation and annihilation. 
On the other hand, if one starts with charged defects of 
infinite energy (which is the case at the point X = 1, 
Y = in the phase diagram in Fig. [2]), then the criti- 
cal dimer state on the bipartite square lattice cannot be 
continuously deformed into an RVB state (see also the 
discussion in Refs. [2^ andlsol). 

There is also an interesting possibility that new phases 
may become possible in our interpolating model, e.g., 
superfluid or supersolid ones. We did not explore those 
possibilities in our analysis, but leave this question for 
future studies. 
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